Accepted by the Astrophysical Journal 

Preprint typeset using L^T^X style cmulateapj v. 10/09/06 



00 

O 
O 

(N 

C 

a 
>—> 

OS 



PHOTOPHORETIC STRUCTURING OF CIRCUMSTELLAR DUST DISKS 
Taku Takeuchi 1 and Oliver Krauss 2 

Accepted by the Astrophysical Journal 

ABSTRACT 

We study dust accumulation by photophoresis in optically thin gas disks. Using formulae of the 
photophoretic force that are applicable for the free molecular regime and for the slip-flow regime, we 
calculate dust accumulation distances as a function of the particle size. It is found that photophoresis 
pushes particles (smaller than 10 cm) outward. For a Sun-like star, these particles are transported to 
0.1 — 100 AU, depending on the particle size, and forms an inner disk. Radiation pressure pushes out 
small particles (< 1 mm) further and forms an extended outer disk. Consequently, an inner hole opens 
inside ~ 0.1 AU. The radius of the inner hole is determined by the condition that the mean free path 
of the gas molecules equals the maximum size of the particles that photophoresis effectively works on 
(100 /im — 10 cm, depending on the dust property). The dust disk structure formed by photophoresis 
can be distinguished from the structure of gas-free dust disk models, because the particle sizes of the 
outer disks are larger, and the inner hole radius depends on the gas density. 

Subject headings: circumstellar matter — planetary systems: formation — solar system: formation 
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1. INTRODUCTION 

Protoplanetary disks are composed of gas and dust. 
At initial stages of their evolution, the gas of protoplan- 
etary disks is as massive as 10~ 3 — 10 _1 Af© (Greaves 

2004) , and small dust particles are mixed with the gas, 
making the disks optically thick at optical wavelengths. 
As the dust particles grow, the number of small particles 
reduces, and at a certain stage, the disks become opti- 
cally thin (Tanaka et al. 2005; Dullcmond & Dominik 

2005) . The amount of gas also decreases, as the gas dis- 
sipates (e.g., Hartmann et al. 1998; Clarke et al. 2001; 
Takeuchi et al. 2005; Alexander et al. 2006a,b). At late 
stages of the disk evolution, the disks become gas-free 
dust disks, which are observed as Vega- type disks. Dur- 
ing the transition from protoplanetary disks to Vega-type 
disks, there may be a stage where the disks are optically 
thin, but their gas component still remains. An example 
is HD 141569A, which is a 5 Myr Herbig Be star (Wein- 
berger et al. 2000) and has an optically thin dust disk 
(Augereau et al. 1999; Weinberger et al. 1999; Fisher 
et al. 2000; Mouillet et al. 2001; Marsh et al. 2002; 
Clampin et al. 2003). The gas component of this system 
has been observed (Zuckerman 1995; Dent et al. 2005; 
Goto et al. 2006) and its amount is estimated to be 
< 60 M ffi (Ardila et al. 2005). 

Dynamics of dust particles in optically thin gas disks is 
of interest in order to investigate the structure of transi- 
tional disks. Krauss & Wurm (2005) considered the mo- 
tion of dust particles in a gas disk that is optically thin 
at optical wavelengths. A dust particle receives stellar 
radiation directly, and the radiation pressure pushes the 
particle outward. In addition to radiation pressure, in- 
teraction between the particle and the surrounding gas 
molecules induces photophoresis, which also pushes the 
particle outward (see also Wurm & Krauss 2006; Krauss 
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et al. 2007). When these outward accelerations act with 
the gas drag on the particle, the particle drifts outward 
(Takeuchi & Artymowicz 2001). In the outer part of the 
disk (> 10 AU), where the mean free path of the gas 
molecules is larger than the particle size (of < 1 m), the 
photophoretic force is proportional to the gas density. 
When a particle moves outward to the point where the 
gas density is no longer high enough to induce a strong 
photophoretic force, the particle's outward motion stops. 
Consequently, the dust particles pile-up at a certain dis- 
tance from the star that is determined by the density pro- 
file of the gas disk. In a gas disk with mass ~ 0.01 M Q , 
particles of 100 ^tm to 10 cm pile-up at several tens of 
AU. This spontaneous ring formation is a characteristic 
feature that is caused by photophoresis in a gaseous dust 
disk. 

Krauss et al. (2007) have demonstrated that pho- 
tophoretic dust motion may be the key process for the 
transition from optically thick protoplanetary disks to 
optically thin circumstellar disks with ring-like dust dis- 
tributions via the stage of transitional disks with an in- 
ner hole and a more or less sharp transition to the outer 
disk that is continuously pushed outward. This outward 
dust migration can also explain the presence of material 
formed close to the sun like chondrules and CAIs in as- 
teroids of the main asteroid belt (Wurm & Krauss 2006) 
or high temperature crystalline silicates in comets from 
the Kuiper belt (Petit et al. 2006; Brownlee et al. 2006; 
Mousis et al. 2007). 

In this paper, we seek to investigate other characteris- 
tic structures formed by photophoresis, especially in the 
inner part of the disk, but at a stage when most of the 
dust has already been built into larger bodies. Krauss 
& Wurm (2005) focus on the dust dynamics in the outer 
part of the disk (> 10 AU) and use a formula for the 
photophoretic force that is relevant for low gas densities. 
In order to investigate the structure of the whole disk, we 
use a photophoresis formula that can be applied to the 
whole range of gas densities. Details of the formula for 
the photophoretic force adopted here are described in fj2l 
Using this formula, we find that photophoresis does not 
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induce dust pile-up at the innermost region of the disk 
(< 0.1 AU), leading to the formation of an inner hole 
(see $3]). At the inner part of the gas disk, the gas den- 
sity may be too high that the gas disk itself is no longer 
optically thin to the stellar radiation, and thus any radi- 
ation effects on the dust such as radiation pressure and 
photophoresis do not work at all. In order to investigate 
whether photophoresis still works at the inner disk, in fJH 
we calculate the optical depth of the disk due to Rayleigh 
scattering of the hydrogen molecules. In M5. It the ther- 
mal relaxation time of dust particles is estimated and 
is compared to the rotation periods induced by gas tur- 
bulence and by photophoresis. In §5.2[ we discuss what 
are the characteristic features that photophoresis forms 
in dust disks, and what are the observable differences 
from the gas-free dust disk structures. In Appendix A, 
a simple calculation deriving the photophoretic force is 
described. 

2. RADIAL DRIFT OF DUST PARTICLES 

2.1. Forces Acting on Dust Particles 

We consider the motion of a dust particle that resides 
in a gas disk around a star. The gas disk is assumed to be 
optically thin in the radial direction. This means that the 
dust particle directly receives the radiation of the central 
star. The stellar radiation that is absorbed or scattered 
by the particle pushes it outward. This outward force, 
which is called radiation pressure (Burns et al. 1979), is 
given by 

Ka 2 Q lad I 



(i) 



where we assume that the particle is spherical and its 
radius is a, Qrad is the efficiency factor of the radiation 
pressure and assumed to be unity in this paper, and c 
is the speed of light. We assume that the particle is on 
the midplane of the gas disk and the distance from the 
central star is r. The flux of the stellar radiation at the 
particle position is 

Aur 2 ' ^ 
where L is the luminosity of the central star. Poynting- 
Robertson (PR) drag can be neglected comparing to gas 
drag. Even in a gas disk whose mass is as small as several 
earth mass, gas drag is much stronger than PR drag as 
discussed below in H3.4I (see also Takeuchi & Artymowicz 
2001). 

If the particle's hemisphere facing the central star be- 
comes hotter than the opposite hemisphere, the particle 
receives the photophoretic force in addition to the radi- 
ation pressure. Because, in general, the particle rotates 
with a certain speed, a considerable temperature gra- 
dient appears only if the thermal relaxation time of the 
particle is much smaller than the rotation period (see dis- 
cussion in §5.1l below). or if the rotation axis is aligned in 
the direction of the light source (Krauss et al. 2007). In 
this paper, we assume that the effect of the particle ro- 
tation can be neglected. In this case, the photophoretic 
force is directed in the radial direction and its value can 
be analytically calculated. In Appendix A, we describe 
estimate of the photophoretic force derived by a simple 
model of cylindrical dust particles, but we use the rigor- 
ous formulae that have been derived in the literature in 
the following discussion. 



We consider a spherical dust particle of radius a. The 
mean temperature Td of the particle is in general differ- 
ent from the surrounding gas temperature T g . In this 
paper, we simply assume that Td — T g . The tempera- 
ture gradient inside the particle is assumed to be small, 
and then the photophoretic force is calculated by solv- 
ing linearized perturbation equation of the temperature. 
In Appendix A, we discuss the cases where the above 
assumptions are not valid. 

If the particle radius, a, is much smaller than the mean 
free path of the surrounding gas molecules, I, i.e., the 
Knudsen number Kn e l/a> 1, then the photophoretic 
force is calculated with the free molecular approximation 
(Hidy & Brock 1967; Mackowski 1989; Beresnev et al. 
1993) and is given by 

F -KaJiPIa 2 

ph ' f " 3[k d T d /a + 4ea SB T$ + aPv T /2] ' 1 > 

where P is the gas pressure, kd is the thermal conduc- 
tivity inside the particle, e is the particle's emissivity of 
thermal radiation and assumed to be unity, o~sb is the 
Stcfan-Boltzmann constant, and vt is the mean thermal 
velocity of the gas molecules. The asymmetry param- 
eter Ji represents how effectively the incident starlight 
induces a temperature gradient inside the particle. If 
the incident light is perfectly absorbed at the particle 
surface, the asymmetry parameter J\ has the maximum 
value 0.5, and it becomes smaller as the temperature 
gradient becomes smaller. The gas molecules that hit 
the particle are assumed to leave the particle with the 
Maxwell velocity distribution of the local surface tem- 
perature of the particle (i.e., the accommodation factor 
a = 1). In many cases, the denominator of equation (J3J) 
is dominated by the first or second terms, i.e., the tem- 
perature gradient in the particle is determined by the 
internal thermal conduction or by the radiative cooling, 
rather than by the thermal conduction to the gas. In this 
case, the photophoretic force is approximated by 

J^Dh f ~ — f= Kn , (4) 

' 3y/2a mo \kd(l + ^ea SB T^a/kd) 

where we used P — n g ksT g and I = 1 / (y2n g a mo \) , and 
n g is the number density of the gas molecules, ks is the 
Boltzmann constant, and cr mo i is the collisional cross sec- 
tion of the gas molecules. In the free molecular regime, 
the photophoretic force is inversely proportional to the 
gas mean free path and proportional then to the gas den- 
sity. 

If the particle size, a, is much larger than the mean free 
path of the gas molecules, I, i.e., the Knudsen number 
Kn <C 1 , then the photophoretic force is calculated with 
the slip-flow approximation (Mackowski 1989). Mack- 
owski (1989) calculated the photophoretic force assum- 
ing that the thermal radiation from the particle surface 
can be neglected. When the thermal radiation is taken 
into account, his formula is modified to (see Appendix B 
for the derivation) 

AV^Cs.hksIa 2 , 
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where 



d = (l+3C m Kn) 
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1(1 + 2C t Kn) 
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and C s , C m , and Ct are the coefficients for the jump 
conditions at the surface and are of order of unity. We 
adopt the values of C s = 1.17, C m = 1.14, and C t = 2.18 
(Mackowski 1989). The thermal conductivity of the gas 
is 

15, 



k n 



—lv T n g k B 



(7) 



where we set the Prandtl number to be 2/3. When d is 
approximated by unity, i.e., Kn <§C 1, Aea^T^a/kd <C 1, 
and kg/kd <C 1, the photophoretic force is proportional 
to the gas mean free path. 

To connect the above two regimes, we introduce the 
following formula 



Kn 2 F ph ,f + F ph ,, 
1 + Kn 2 



(8) 



In the limit of Kn <C 1, F p h approaches i^ph.s, and in the 
other limit (Kn S> 1), it approaches f p h,f- 

If the particle rotation cannot be neglected, there are 
two effects we must consider. First, the temperature gra- 
dient in the particle becomes smaller as the particle ro- 
tates faster. This effectively reduces the asymmetry fac- 
tor J\. In the limit of rapid rotation, the asymmetry 
factor J\ approaches zero. Second, the temperature gra- 
dient generally does not coincide with the direction of the 
light source, i.e., the radial direction. The photophoretic 
force has a component in the azimuthal direction, and 
can give or take away angular momentum of the parti- 
cle, depending on the orientation of the rotation. The 
particle drifts outward or inward, and the direction of 
the drift reverses when the orientation of the rotation 
changes. This is similar to the Yarkovsky effect acting 
on asteroids (Bottke et al. 2002), although, in the case 
of photophoresis, the agents of momentum transfer are 
the gas molecules and not emitted photons. The particle 
radial drift due to this effect is stochastic, if the variation 
in the rotational orientation occurs frequently before the 
particle travels a large distance in the radial direction. 
Though the study of photophoretic Yarkovsky effect is 
important, we leave this for future investigations, and 
in this paper we assume that the photophoretic force is 
always directed to the radial direction. 

For further discussion, it is convenient to normalize the 
radiative forces given above by the gravitational force 



F G , d 



GMm d 



(9) 



where G is the gravitational constant, M is the mass of 
the central star, and m d is the mass of the particle. We 
introduce a normalization factor 



= A-ad + Pph , 



where 



and 



/3rad 



-^rad 



F G ,d 
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In the free molecular regime (Kn 3> 1), 
aJxksL 

Pph 



16V2irGMa mol k d p d l{l + 4ea S BT$a/k d ) 



(10) 
(11) 
(12) 

(13) 



At large distances from the star, the dust temperature is 
so low that ieasBTdd/kd <C 1, and that f3 p h is indepen- 
dent of the particle radius a. 

Next we consider the forces acting on the gas. The 
gravity acting on a unit volume of the gas is 



F G , g 



GMpg 



(14) 



where p g is the mass density of the gas. The pressure 
gradient force on a unit volume is 



F^p — — 



dP 
dr 



(15) 



We introduce a normalization factor rj in a similar way 
to the dust particle, 



V 



Fvi 



F G,g 



(16) 



2.2. Gas Disk Models and Dust Parameters 



For simplicity, we assume that the gas disk has a 
power- law profile of the temperature, T g , in the radial 
direction r, and is isothermal in the vertical direction z, 
i.e., it is written as 

T g {r)=T gfi r-« 



AU ' 



(17) 

where the subscript "0" denotes quantities at 1 AU, and 
the non-dimensional quantity tau is the radius in AU. 
The gas density, p g , is also assumed to have a power-law 
profile in r and be in a hydrostatic equilibrium in z, i.e., 



(18) 



p g {r, z) = Pgfirjp exp - — 



2h* 



where the disk scale height is V2h g . The isothermal 
sound speed of the disk is c s = c r A ^ 2 , and h g is written 

by 



Mr) = = AorfcT*)' 3 



(19) 



where Ok = (GM/r 3 ) 1 ' 2 is the Keplerian angular fre- 
quency at the disk midplane. The gas surface density 



is 



S fl (r) 



+ 00 



p g dz = V2Trp h r A ^ , 



(20) 



where p = p m + (q— 3)/2. We adopt the following fiducial 
parameters for the central star and the gas temperature 
profile: M = 1 M Q , L = 1 L , T gS) = 278 K, h = 3.33 x 
10~ 2 AU, and q = \. The mean mass of the gas molecules 
is 2.34t7j h , where m H is the mass of a hydrogen atom 
(Nakagawa et al. 1986). The collisional cross section of 
the gas molecules is cr mo ] ~ 2 x 10 -15 cm 2 (p. 228 in 
Chapman & Cowling 1970). 

We consider disks in which most of small dust grains 
have been removed and consequently the disks have be- 
come optically thin even in the radial direction. The 
dust opacity in optical wavelength arises mainly from 
particles smaller than 10 pm (see Fig. 4 of Miyake & 
Nakagawa 1993). In order for the disk to be optically 
thin, the column density of small grains (< 10 pm) from 
the star must be smaller than 10 -3 g cm~ 2 . Removal 
of dust grains probably occurs through coagulation of 
the grains into planetesimals (Weidenschilling & Cuzzi 
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1993; Tanaka et al. 2005; Dullemond & Dominik 2005; 
Nomura & Nakagawa 2006). In such disks, the amount 
of the gas component has likely also reduced. Actually 
most Vega-type stars do not have detectable gas com- 
ponents (Liseau & Artymowicz 1998; Greaves et al 2000; 
Coulson et al. 2004; Chen & Kamp 2004), and only a few 
objects have gas as massive as several tens earth masses. 
For example, HD 141569A has a gas mass of < 60M e 
(Zuckerman et al. 1995; Ardila et al. 2005). Therefore, 
in model A, we consider a disk in which the amount of 
the gas inside 100 AU is M g = 2.4 x 10" 5 Af Q = 7.9M e . 
(The model disk extends over 100 AU. We specify the 
disk mass as the mass inside 100 AU.) We take a sur- 
face density profile proposed by Hayashi et al. (1985) in 
which the power-law index is p = 1.5. The gas surface 
density is S g = 1.7 r^y g cm -2 . In addition to model 
A, we also consider a more massive disk to investigate a 
case in which the disk opacity due to the dust has been 
significantly reduced before considerable gas dissipation 
occurs. In model B, the gas amount inside 100 AU is 
M g = 1.2 x 10~ 2 M Q . Further, we assume the density 
profile is more gradual than model A and p = 0.5. This 
gentle slope of the density profile assures that the gas 
density in the innermost part of the disk (say ~ 0.1 AU) 
does not become so high that the Rayleigh scattering by 
molecular hydrogen makes the gas disk itself become op- 
tically thick to the starlight (see §4.11 below) . The gas 
surface density in model B is T, g = 25 r^u' 5 g cm" 2 . 

For dust particles, we use the following fidu- 
cial parameters: The thermal conductivity kd = 
10 2 erg s" 1 cm" 1 K" 1 , the particle bulk density pd — 
1 g cm -3 , and the asymmetry factor J\ — 0.5, the emis- 
sivity £ = 1, and perfect accommodation of gas molecules 
a = 1. The adopted value of the thermal conductivity 
is typical for porous aggregates (Presley & Christensen 
1997). If the dust particles are more compact, the ther- 
mal conductivity can be much higher. 

2.3. Radial Motion of Dust Particles 

Both the radiation pressure and the photophoretic 
force direct outward. If the sum of these forces is stronger 
than the gravity of the central star, i.e., /3 > 1, then the 
particle is not bound to the central star and is ejected 
from the system. If (3 < 1, the particle is bound to the 
central star, but it feels a net force weaker than the grav- 
ity of the central star, and its orbital motion is slower 
than Keplerian. Assuming the particle's orbit is circular 
and neglecting any contribution of the gas drag force, the 
azimuthal velocity of the particle is 

«M = (1-/3) 1/2 «K , (21) 

where «k = (GM / r) 1 / 2 is the Keplerian velocity. The az- 
imuthal velocity of the gas is also sub-Keplerian because 
the pressure gradient force usually directs outward, and 
is given by 

ve, g = (l-v) 1/2 VK ■ (22) 

Now we consider the effect of the gas drag force act- 
ing on a particle that is bound to the central star (i.e., 
(3 < 1). If the particle orbits faster than the surrounding 
gas (vg t d > vg t g), then the gas drag transfers the angu- 
lar momentum of the particle to the surrounding gas. 
Consequently, the particle loses angular momentum and 
drifts radially to the central star. From equations (|21[) 



and (|22|) , it is seen that this inward drift occurs if f3 < rj. 
On the other hand, if the particle orbits slower than the 
surrounding gas {vg t d < ve,g, i.e., (3 > rj), then the par- 
ticle gets angular momentum from the surrounding gas, 
and drifts radially away from the central star. When 
the particle happens to be on the orbit in which its az- 
imuthal velocity equals the gas azimuthal velocity, there 
is no angular momentum transfer by gas drag, and the 
particle stays on the equilibrium orbit. More detailed 
discussions including the derivation of the drift velocity 
were written in Takeuchi & Artymowicz (2001). From 
their equation (26) for the drift velocity, neglecting the 
term of the Poynting- Robertson drag, (3 C T S , gives 



Vr,d = 



T 8 + T s ~ 



-v K 



where the non-dimensional stopping time is 



p g rv T 



(23) 



(24) 



[the factor 4/3 difference between eq. (10) of Takeuchi & 
Artymowicz (2001) and eq. ([24]) comes from the different 
definition of i>t, see their eq. (6)] 

We can know the direction of the particle's radial drift 
by seeing if (3 is larger than rj or not. Figure [TJi shows 
the ratio (3/rj against the distance from the central star 
for model A. If a particle has (3 larger than r\ of the sur- 
rounding gas, it drifts away from the star, and vice versa. 
The ratio of radiation pressure to gravity, (3 Ta d, which is 
inversely proportional to the particle radius a, is con- 
stant with the distance from the star r (see eq. [H]). On 
the other hand, as long as the gas mean free path I is 
larger than the particle radius a (free molecular regime), 
the photophoretic force is inversely proportional to the 
Knudscn number Kn (see eq. [3]), and hence decreases 
with r. Thus, at large distances from the central star, 
the radiation pressure dominates the photophoretic force, 
i.e., (3 ~ /3 ra dj while at small distances, the photophoretic 
force dominates and (3 ~ (3 p h- The transition occurs at 1 
AU for 10 /iin particles, 3 AU for 100 pm particles, and 7 
AU for 1 mm particles. The curves in Figure [TJz. bend at 
these distances and the slopes become gentler at larger 
distances. The slopes of the curves, d((3 / 'rj) / dr , are deter- 
mined only by drj/dr at large distances, because /3 rac i is 
constant with r, and have the same value for differently 
sized particles. Inside the transition distances, where the 
photophoretic force dominates, /3(« /3 P h) increases with 
decreasing the Knudsen number Kn, or decreasing the 
distance from the central star. This increase in (3 con- 
tinues until the mean free path of the gas molecules I 
becomes smaller than the particle radius a, i.e., until Kn 
becomes less than unity. At the innermost part of the 
disk in which Kn < 1 (slip- flow regime photophoresis), (3 
is proportional to Kn (see eq. [5]) and decreases with de- 
creasing r. The transition from the free molecular regime 
photophoresis to the slip-flow regime photophoresis oc- 
curs at the distance where Kn = 1, and at this transition 
distance, (3/rj has the maximum value. For 100 pm par- 
ticles, (3/rj has a peak at 0.02 AU, and for larger particles 
the peak position shifts outward. It is 0.2 AU for 10 cm 
particles. 

It is seen from Figure [Th that particles smaller than 
10 cm have a region where (3 is greater than 77. In such 
a region the particles drift outward. The outward drift 
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Fig. 1. — Ratio of the normalized outward force on the dust 
particle, f3, with regard to that on the gas, rj. The particle size 
is varied from 10 /mi to 1 m. The dotted line represents /3 = 
r\. The particles drift outward if f3/r] > 1, and drift inward if 
f3/r] < 1. (a) For the gas disk model A (S s = 1.7 r^J' 5 g cm' 2 , 
M g = 2.4 x lCr 5 M Q = 7.9M® inside 100 AU). (6) Model B (E 9 = 
25 r~°- 5 g cm" 2 , M g = 1.2 X 10- 2 M Q ). 

of the particles stops at the locations where (3 = r\. At 
these locations, the particles that drift inward from the 
outer region also stop. Hence, the particles accumulate 
at the location where (3 = r\ and d((3/n)/dr < 0, leading 
to the formation of a dust ring as proposed by Krauss & 
Wurm (2005). Particle accumulation does not occur at 
the location where (3 = r\ but d(f3/r])/dr > 0, because the 
particles drift away from such locations. Note that, at 
the locations where particles accumulate, the Knudsen 
number Kn is always larger than unity, and thus the 
photophoretic force is expressed in its form for the free 
molecular regime. The location of particle accumulation 
is at 5 AU for 100 /im particles, and moves closer to the 
central star as the particle radius increases. Particles of 
10 cm accumulate at 0.2 AU. If the particles are larger 
than 10 cm, they always drift toward the central star 
and cannot accumulate anywhere. In summary, particle 
accumulation occurs if fl/rj > 1 somewhere in the disk. 
The accumulation location is specified by the conditions 
(3 = r], d(0/r])/dr < 0, and Kn > 1. 

For 10 cm and 1 m particles, the curves in Figure \^p, 
again bend at 0.02 — 0.05 AU. Inside these distances, 
the radiation pressure becomes stronger than the pho- 



tophoretic force because the Knudsen number is too low. 
Thus, [3 keeps constant with r in this region, while r\ 
increases, resulting in a decrease in (3/n. 

A similar result is obtained for model B (Fig. [IJ)). 
In this model, particles of all sizes are driven farther 
outward compared to model A. This is because the gas 
density is higher in the model B disk, resulting in an 
increase in /3 p h- In model B, it is also apparent that 
100 /im — 1 cm particles accumulate in a narrow region 
of 20 — 30 AU. This pile-up of particles in a narrow ring 
has been pointed out by Krauss & Wurm (2005). It is 
the consequence of /3 p h being independent of the particle 
radius at large distances (eq. [32]). 

Note that accumulation of outward migrating particles 
to the ring radius do not increase the optical depth to 
the star. Formation of the dust ring through clearing 
the inner dust disk does not suppress further evolution 
by photophoresis. However, after clearing the inner dust, 
the optical depth of the ring increases as particles from 
the outer disk accumulate. When the dust ring becomes 
optically thick, photophoresis is weakened and then the 
ring shrinks. In this paper, we consider the stage where 
the dust ring is optically thin. 

3. PARTICLE ACCUMULATION RADIUS 

As discussed in the above section, particles accumu- 
late at the equilibrium locations where (3 = rj and 
d(f3/n)/dr < 0. In Figure [2J3, and[2f) the equilibrium dis- 
tances for model A and model B, respectively, are plot- 
ted as functions of the particle radius. The line labeled 
"fiducial" represents the equilibrium distances calculated 
with the parameters written in ^2.21 The equilibrium dis- 
tances for other parameters are also plotted. We vary the 
value of one parameter with keeping the other parame- 
ters being the same as the fiducial model. Higher ther- 
mal conductivities kd = 10 5 and 10 7 erg s _1 cm -1 K _1 , 
a lower bulk density pd — 0.1 g cm~ 3 , a lower efficiency 
of photophoresis Ji = 0.5 x 10~ 2 , and a higher luminos- 
ity of the central star L — 20 Lq are investigated. We 
also calculate the case in which photophoresis does not 
work at all (Ji = 0). 

3.1. Dust Concentration and Size Segregation 

At the accumulation locations, photophoresis is un- 
der the free molecular regime (Kn > 1), and as seen 
from equation (fT3")) . /3 ph converges to a certain value in 
the limit of small a. On the other hand, /3 ra d is in- 
versely proportional to a (neglecting the size dependence 
of Qrad), and thus radiation pressure is more efficient 
than photophoresis for smaller particles. In the condi- 
tion (3 = r\ that determines the accumulation location, 
(3 is controlled by radiation pressure for small particles, 
and by photophoresis for large particles. The particle 
size where the transition from one behavior to the other 
occurs depends on the parameters that are relevant for 
the strength of the photophoretic force, i.e., kd and Ji. 
This can be seen in Figure [2)2 where the curves for differ- 
ent parameters depart from the straight line for J\ = at 
different particle sizes. For the fiducial dust parameters 
in model A (the solid line in Fig. Wp), the accumulation 
locations of particles larger than 100 /jm are controlled 
by photophoresis, while those of smaller particles are con- 
trolled by radiation pressure. 

Particles that are subject to photophoresis have the 
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Fig. 2. — Equilibrium distances, where j3 = t], for the 
dust particles of the various physical properties: fiducial pa- 
rameters (krf = 10 2 erg s _1 cm -1 K" 1 , p& = 1 g cm -3 , and 
Ji = 0.5; solid line), higher thermal conductivity (kj = 10 5 and 
10 7 erg s" 1 cm" 1 K — 1 ; short dashed lines), smaller efficiency of 
photophoresis (J\ = 0.5 X 10 — 2 and 0; long dashed lines), lower 
bulk density (p^ = 0.1 g cm" 3 ; dot-dashed line), and more lu- 
minous star L = 2OZ/0 (three dots-dashed line). The dotted line 
shows the distances where the Knudscn number is unity. Above 
(below) the dotted line, Kn > 1 (Kn < 1). (a) Model A. (6) Model 
B. 

tendency to accumulate in a narrow ring region, provided 
that all particles have the same physical properties, i.e., 
they have the same values of pd, kd, and J\. For example, 
in the fiducial model (the solid line in Fig. [2b), particles 
of 200 /im — 2 mm accumulate in a region of 2 — 4 AU, 
and 300 /im — 10 cm particles of the higher conductivity 
k d = 10 5 erg s -1 cm" 1 K" 1 accumulate at 0.2 — 0.5 AU. 
Hence, a significant influence of photophoresis is to con- 
dense the particles of a certain size interval. For small 
particles, this size interval is limited by the increasing 
influence of radiation pressure. For large particles, the 
limit is determined by the increasing importance of ther- 
mal emission and conduction to the gas for the thermal 
relaxation and the approach to Kn = 1. The location 
of accumulation decreases with increasing thermal con- 
ductivity. Thus, particles with different thermal proper- 
ties can be separated in the disk by this process. While 
kd = 10 2 erg s _1 cm" 1 K" 1 corresponds to very porous 
dust aggregates, values around 10 5 erg s _1 cm" 1 K _1 
are appropriate for compact rocky or glassy grains and 



10 7 erg s _1 cm" 1 K" 1 for purely metallic particles. As 
discussed by Wurm & Krauss (2006), compact particles 
like chondrules or CAIs might have been accumulated in 
the region of the asteroid belt whereas porous dust aggre- 
gates were driven farther outwards by photophoresis in 
the forming solar system (see Fig. [2j) for the disk model 
as massive as the solar nebula). 

In actual disks, the properties of particles are probably 
not homogeneous, but the dust disks are composed of a 
wide variety of particles. Even particles of the same size 
therefore spread over a certain range of distances, and 
the dust ring formed by photophoresis has some extent 
in the radial direction. We still expect, however, particle 
concentration in the disk happens to some degree. Even 
if the parameter range of the particle properties is as wide 
as that considered in Figure [2bj we clearly see two pop- 
ulations of dust particles. All the particles at distances 
larger than 10 AU are smaller than 1 mm and do not 
show any sign of concentration. The dust population in- 
side 10 AU is composed of particles larger than 100 /im, 
and is confined in a region between 0.1 — 10 AU. We call 
this region as an inner disk. These two populations are 
divided by the condition which of the radiation pressure 
or the photophoretic force is larger than the other. The 
population of small particles is subject to radiation pres- 
sure, and the other is subject to photophoresis. Hence, 
the effect of photophoresis is to discriminate large parti- 
cles from small particles and to confine them to the inner 
region of the disk. 

3.2. Inner Holes 

Another significant feature of Figure [2b is that in the 
region below the dotted line, which represents the area 
where Kn < 1, there is no equilibrium location where 
particle accumulation through photophoresis occurs. As 
discussed in £12.3i particle accumulation by photophoresis 
always occurs at locations where Kn > 1. If photophore- 
sis does not work (Ji=0), it is seen in FigureGb that par- 
ticles can accumulate even in the region where Kn < 1. 
This feature will be a possible observational evidence for 
showing that photophoresis actually occurs in gas disks, 
though it requires observations of the innermost part of 
the disk (r < 0.1 AU) with a very high spatial resolution. 
If photophoresis actually works, we would find that an 
inner hole opens in the dust disk, and that its radius is de- 
termined by the condition Kn = 1. If photophoresis does 
not work, on the other hand, the dust particles would ex- 
ist even in a region where Kn < 1. An actual dust disk 
is a mixture of particles of various physical properties. If 
the dust particles have the range of properties that are 
considered in Figure [21 the inner radius (that is deter- 
mined by the particles of kd = 10 7 erg s^ 1 cm" 1 K" 1 or 
Ji = 0.5 x 10 -2 ) is ~ 0.1 AU. This inner radius may be 
as small as 0.03 AU, if there are considerable numbers of 
particles on that photophoresis works very weakly (i.e., 
particles of very small J± , but in that case the inner disk 
clearing is probably suppressed as shown by the J\ = 
line.) 

We stress that the Knudsen number Kn is only a func- 
tion of the dust particle size a and the mean free path of 
the gas molecules Thus, the inner radii of dust disks do 
not directly reflect the dust particles' physical properties 
such as the thermal conductivity kd, the bulk density pd, 
and the efficiency of photophoresis J\ . (The particle size 
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Fig. 3. — Timescales, T m j g , for the dust particles to migrate to 
the accumulation locations. The dust parameters are the fiducial 
values. 

at the inner edge does depend on the particle proper- 
ties.) The fact that the inner radius of the dust disk is 
determined only by the particle size considerably helps 
us when we observationally test the idea of inner hole 
opening by photophorcsis. On the other hand, it is of 
great importance to have an independent measurement 
of the local gas density to decide whether the condition 
Kn = 1 is fulfilled, which might still be a challenging task 
for modern observational techniques. 

3.3. Massive Disks 

Figure [2> shows the equilibrium distances with the 
same parameter combinations as in Figure [5Jz but for the 
disk model B. In model B, where the disk mass (inside 
100 AU) is 500 times larger than that of model A, 1 mm 
to 10 cm particles may extend farther from the central 
star than in model A due to photophoresis (up to ~ 100 
AU). The inner radius of the dust disk locates at ~ 0.1 
AU (for the particles of kd — 10 7 erg s^ 1 cm -1 K _1 or 
Ji = 0.5 x 10 -2 ), and is similar to that in model A. This 
is just because we take such a density distribution for 
model B that the distances where Kn equals unity for 
the particle size a = 1 mm — 1 m do not considerably 
differ from those in model A. This choice of the den- 
sity distribution in model B comes from the requirement 
that the gas disk must be optically thin and we take a 
marginally optically thin disk (see H4.ll below) . 

3.4. Timescale of Structure Formation 

The timescale for the dust particles to migrate to the 
accumulation locations is r m i g = \r cq /v r ,d\, where r oq is 
the equilibrium distance shown in Figure [5] and v r ,d is 
the radial drift velocity. In equation (|23[) . (3 ~ rj near the 
equilibrium distance r eq , and thus the drift velocity v r ^d 
is estimated as u r ,d ~ f] v K/(T s + T~ ). The migration 
timescale is then r m i g ~ (T s + Tr 1 ) / '(ryf2x)) where the 
value is calculated at r oq . Figure [3] shows the migration 
timescales for the fiducial dust parameters of model A 
and B. For both models, the timescales are shorter than 
10 6 — 10 7 yr, which is the expected age of the transitional 
disks, except for small particles (a < 100 /im) in model B. 
Hence, the disk structure and the inner hole are expected 
to form quickly when the disks have become transparent 
to the starlight. (In disks more tenuous than model A, 



photophoresis effectively works only at the dense inner- 
most part of the disk, and thus the dust accumulation 
location approaches the star. For example, in a disk 
10~ 3 times more tenuous than model A, i.e., in a disk 
of M g ~ 1O _2 M0, particle accumulation occurs only in- 
side 0.1 AU. The migration timescale for 100 /urn — 1 cm 
particles at < 0.1 AU is still less than 10 4 yr.) 

The migration timescale in Figure [3] is much shorter 
than the orbital decay time due to Pointing-Robertson 
(PR) drag. For example, 100 /jm particles in model A 
accumulate at 4 AU in r m i g = 2 x 10 3 yr. The particles 
have /? ra( j = 6x 10~ 3 and the timescale due to PR drag 
is r PR = 400(r/l AU) 2 /?"^ yr = 10 6 yr (Burns et al. 
1979), and much longer than r m j g . Both in model A and 
B, it is shown that T m ; g -C tpr for any particle size. This 
means that PR drag is much weaker than gas drag and 
can be neglected. 

4. OPACITY OF THE DISK GAS 

In the previous sections, we assumed that the disk is 
optically thin even in the radial direction and dust parti- 
cles directly receive the light from the central star. This 
assumption requires that both the dust and gas opacities 
of the disk must be small enough. We consider a stage in 
which most of the dust grains are packed together into 
planetesimals or larger objects and the disk opacity due 
to the dust has become negligible. Such a situation has 
probably been realized in some Vega- type stars. The gas 
disk also must be transparent, though photophoresis re- 
quires a sufficient amount of gas. In this subsection, we 
discuss whether the gas disk that is dense enough for 
photophoresis to work is transparent to the starlight. 

We assume that the gas disk is mostly composed of 
hydrogen molecules. Thus, Rayleigh scattering by H2 
is a possible and probably dominant mechanism of light 
extinction in the disk. In this subsection, we consider 
the gas opacity due to Rayleigh scattering by H 2 and 
estimate the maximum gas density at which the gas disk 
is marginally opaque to the light of the central star. Our 
estimate gives the lower limit of the gas opacity, because 
only Rayleigh scattering by H2 is taken into account, 
and thus the derived maximum gas density should be 
considered as an upper limit. A transparent gas disk 
in which photophoresis works must be less dense than 
our estimate. If other extinction mechanisms work more 
effectively than the Rayleigh scattering by H2 , this upper 
limit of the density becomes smaller. In ij4.2l below, we 
briefly discuss other extinction mechanisms. 

4.1. Rayleigh Scattering by Hydrogen Molecules 

The cross section of Rayleigh scattering by H2 per 
molecule is 

cr H2 = 8.4909 x 10~ 29 ( rr~— — ^ cm 2 , (25) 
\ 1 /im J 

where A is the wavelength of the light (Tsuji 1966). In 
this subsection, we make an order of magnitude argu- 
ment, and simply estimate the opacity around the central 
wavelength of the central star's spectrum, A = 0.5 ^m. 
The stellar light source is assumed to be point sources, 
which are located slightly above and beneath the origin, 
= ±0.424i?„, where i?* is the radius of the central star 
(Hollenbach et al. 1994) and we assume i?* = 1.0 Rq. 



Takeuchi & Krauss 



1.0 










depth 


\w ^ Mode 


B 




o 
l> 








<x 
O 


/ Mode 


A 




0.1 









0.01 0.10 1.00 10.00 

Distance from the star (AU) 



100.00 




0.01 0.10 1.00 10.00 100.00 

Distance from the star (AU) 



Fig. 4. — Optical depth th, due to Rayleigh scattering of H2 
from the typical point of the star (0, 0.424i? t ) to the midplane of 
the disk at (r, 0). The dotted line represents tjj, = 1. 

The optical depth from the star, (0, z*), to a dust particle 
on the midplane, (r, 0), is 

th 2 = J <7H 2 n g ds , (26) 

where n g is the number density of the gas molecules (we 
assume all the gas is composed of H2 for simplicity) . 

Figure [4] shows the optical depths from the central star 
as functions of the distance r. In both models A and B, 
the disks are optically thin (except for r < 0.03 AU in 
model A). Note that the disk mass (inside 100 AU) of 
model B is 500 times larger than that of model A, but 
the optical depths from the star are not much different. 
This is because the opacities of the disks are mainly de- 
termined by the density at ~ 0.1 AU and the densities 
at 0.1 AU in model A and B do not differ much (Fig. 
[5j. The starlight that illuminates a particle located at 
a distance r > 0.1 AU gets into the disk at a distance 
~ 0.1 AU, because the disk scale height at 0.1 AU is 
y/2h g « 2 x 10" 3 AU w 0.4R* w z*. Model B has the 
density profile S s cx r -0 5 and most of the gas mass is 
distributed in the outermost part of the disk, while the 
density profile of model A is more steep (E s cx r~ 15 ) 
and the disk mass is more concentrated in the innermost 
region. Thus, even though the disk mass of model B is 
500 times larger than that of model A, the densities at 
0.1 AU in model A and B are similar. For particles lo- 
cated much farther than 0.1 AU, the optical depth from 
the central star does not change with the distance r. 

For particles inside ~ 0.1 AU, the location where the 
starlight enters the disk moves toward the central star 
as the particle position gets closer to the star. In model 
A, the optical depth is determined by the gas density 
at where the starlight enters the disk, which rapidly in- 
creases with decreasing the distance r. Consequently, 
the optical depth increases with decreasing the particle 
distance. In model B, the concentration of the gas at 
the innermost part of the disk is moderate. Thus, the 
decrease in the path length of the starlight (which is ap- 
proximately proportional to r) dominates the increase in 
the gas density at where the starlight enters the disk. 
The optical depth from the star decreases with decreas- 
ing the particle distance in model B. 

We estimate the maximum gas density of a disk for it 



Fig. 5. — Comparison of the density profiles of models A and B 
(dashed lines) to the maximum gas density required to be optically 
thin to the central star (solid line). The dotted lines show the 
density profiles that have the same amount of mass as model B 
(inside 100 AU) but have different slopes (S 9 oc r~ 5 and S 9 oc 
r -1.0), 

to be optically thin. The optical depth to the star from 
(r, 0) at the disk midplane is estimated as t w an 2 n g s, 
where the number density of hydrogen molecules n g is es- 
timated at (r, 0), and s is the path length of the starlight. 
The path length s is the distance between (r, 0) and 
the point where the starlight enters the disk. The lo- 
cation of the latter point is approximately calculated as 
the point where the line connecting (0,0) and (r, \f2h g ) 
crosses the line connecting (0,2*) and (r, 0). The condi- 
tion r < 1 reduces to the condition for the gas surface 
density E 9 < \ // 2TTh g mn 2 /(aH 2 s), where mn 2 is the mass 
of a hydrogen molecule. The solid line in Figure [5] shows 
the upper limit of the gas surface density for the optical 
depth to the star to be less than unity. We see that for 
r > 0.02 AU both the surface densities of model A and B 
(dashed lines) are less than the critical density for r = 1 
(solid line). For comparison, the surface densities of the 
disks that have the same amount of mass as model B 
(inside 100 AU) but have different slopes of the density 
distribution (£ g cx r~ 15 and E ff cx r~ 10 , dotted lines). 
It is seen that, if the disk is centrally concentrated, the 
disk is optically thick at the inner region of the disk. 

The condition t < 1 reduces also to the condition for 
the mean free path of the gas molecules. The mean free 
path of the gas must be larger than the value shown as 
the solid line in Figure [5] for the disk to be optically thin. 
This minimum value is also interpreted as the minimum 
size of dust particles that compose the inner edge of the 
dust disk, because the location of the inner edge is de- 
termined by the condition Kn = 1. Thus, in order to 
observationally test whether the inner hole is produced 
by photophoresis, we have to detect particle accumula- 
tion of sizes larger than the value indicated in Figure 
[6l Suppose that, for example, a dust disk has an inner 
hole of radius 0.04 AU. If this inner hole is formed by 
photophoresis, the gas disk at 0.04 AU must be optically 
thin, and from Figure O the mean free path at 0.04 AU 
must be larger than 1 mm. Because only particles larger 
than the mean free path can accumulate through pho- 
tophoresis, the particles at the inner edge must be larger 
than 1 mm. This minimum particle size becomes larger 
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Fig. 6. — Mean free path of the hydrogen molecules in the 
critical disk where the optical depth to the central star is unity. 

as the inner radius of the dust disk becomes larger. If 
the inner radius is 1 AU, the particle size must be larger 
than 10 cm. 

4.2. Other Opacity Sources 

1. Bound-free and free- free absorption of negative hy- 
drogen ions — The cross section of bound-free and free- 
free absorption is proportional to the electron pressure. 
Thus, this absorption works in environments where the 
ionization degree is high, and is efficient only at the in- 
nermost region of the disk, provided that the disk gas 
is thermally ionized. In both models A and B, the tem- 
perature at 0.02 AU is T w 2000 K, and the gas pres- 
sure is Ph 2 ~ 10 3 dyn cm -2 (the molecular hydrogens 
dominate the gas pressure). The partial pressures of 
electrons and of hydrogen atoms are P e /Pn 2 ~ 10~ 7 
and Ph/Ph 2 ~ 3 x 10~ 2 , respectively (Fig. 1 of Tsuji 
1966; the solar abundance is assumed). At the op- 
tical wavelength A = 0.5 /im, the bound-free absorp- 
tion dominates the free-free absorption, and the absorp- 
tion cross section per hydrogen molecule is 0h- ,&/+// ~ 
10~ 23 P e P H /fH 2 cm 4 dyn- 1 - 3 x 10~ 29 cm 2 , where we 
used the value of Table 2 of Tsuji (1966). Since the cross 
section of the bound-free absorption sensitively depends 
on A and T, we took the highest value around A rs 0.5 /im 
and T ss 2000 K in the Tsuji's Table as a representative 
value (we take the value for A = 0.5 /im and T = 1938 
K). From equation (|25p. the cross section of Rayleigh 
scattering for A = 0.5 /im is <th 2 = 1-36 x 10~ 27 cm 2 . 
The cross section of bound-free and free-free absorption 
is much smaller than that of Rayleigh scattering even at 
the innermost part (r = 0.02 AU) of disk model A and 
B, and thus can be ignored. (Recall that stellar photons 
penetrate the disk at ~ 0.1 AU for calculating the optical 
depth tor» 0.1 AU.) 

2. Free-free absorption of negative hydrogen molecules 
- This absorption works only at the innermost part 

of the disks. Thus, we estimate the cross section 
at 0.02 AU of the model disks, which is <j h - ^ ~ 

10~ 27 P e cm 4 dyn -1 ~ 10~ 31 cm 2 , where we used equa- 
tion (8) of Tsuji (1966). This value is much smaller than 
that of Rayleigh scattering, and free-free absorption of 
negative hydrogen molecules can be ignored. 

3. Collision-induced absorption of hydrogen molecules 



- This absorption works efficiently if the gas pres- 
sure is high. The cross section is expressed as 
ctcia — kciAn-g/n^ where kciA is the opacity in the 
unit cm -1 amagat -2 and = 2.69 x 10 19 cm -3 is 
Loschmidt's number. We consider the innermost part of 
the disk. The gas number density at 0.02 AU of model 
A and B is n g ~ 10 16 cm -3 . At the optical wavelength 
A = 0.5 /xm, the opacity fcciA due to H2 — H2 and H2 — He 
collisions is less than 10 -11 cm -1 amagat -2 as long as 
T < 4000 K (Fig. 1 and 2 of Borysow et al. 1997). Thus, 
ccia is smaller than 10" cm 2 , and can be ignored in 
comparison to Rayleigh scattering. 

4. Thomson scattering — The cross section of Thom- 
son scattering per hydrogen molecule is a e — 6.65 x 
10 _25 P e /PH 2 cm 2 . This can be ignored in comparison 
to Rayleigh scattering, as long as P e /Pn 2 <C 10~ 3 . 

5. Molecular absorption bands — At the innermost 
part of the disk, where the temperature is high enough to 
vaporize refractory elements from the dust, gas molecules 
such as TiO and VO contribute to the gas opacity at 
the optical wavelengths (Tsuji 1971; Alexander & Fergu- 
son 1994; Ferguson et al. 2005). The rotation- vibration 
lines of such molecules overlap and make a band struc- 
ture veiling the optical wavelengths. We calculated the 
gas opacity of model disks A and B due to the molec- 
ular absorption bands. In the calculation, we used the 
code developed by Tsuji (2002) for calculating the opac- 
ity of cool stellar atmospheres, and assumed the solar 
abundance of the gas and dust mixture. At the optical 
wavelengths (A ~ 0.5 /im), absorption by VO becomes 
stronger than Rayleigh scattering of H2 for T > 1300 K, 
and absorption by TiO becomes stronger for T > 1400 
K. Thus, in the disk inside 0.05 AU of models A and B, 
where the temperature is higher than 1300 K, the gas 
opacity is dominated by the molecular absorption bands. 
Our estimate of the gas opacity in £14.11 is not appro- 
priate for the gas inside 0.05 AU. For dust particles at 
r "3> 0.1 AU, the rays from the central star enter the 
disk at ~ 0.1 AU, where the gas opacity is dominated by 
Rayleigh scattering, and contribution of the molecular 
absorption bands can be neglected. Most of refractory 
elements may be confined in large bodies such as plan- 
etesimals and are removed from the gas phase. In such 
cases, contribution of the molecular absorption bands is 
reduced. 

5. DISCUSSION AND SUMMARY 

5.1. Thermal Relaxation and Rotation Times of Dust 
Particles 

As mentioned in £12.11 if particles rotate rapidly, pho- 
tophoresis does not work. Even if particle rotation is 
not considerably rapid compared to the thermal relax- 
ation time of the particle, photophoretic Yarkovsky ef- 
fect may prevent the particle accumulation. Such cases 
where photophoresis does not work are represented by 
the model of J\ = 0, where the incident starlight does not 
cause any temperature gradient in the particles. If pho- 
tophoresis is suppressed by any reason, particles can re- 
side in the region where Kn < 1 (Fig. [2]), and such cases 
can be observationally distinguished from the disks in 
which photophoresis clears the dust in the innermost re- 
gion. Hence, whether particles in gas disks rotate rapidly 
or not can be observationally investigated by checking 
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Fig. 7. — Thermal relaxation times of cylindrical particles of 
various sizes are plotted by the solid lines against the distance 
from the star. The dashed line shows the turnover time, r e( j, of the 
smallest eddies for atur = 10 -2 . 

whether the dust in the region of Kn < 1 is cleared or 
not. 

In this subsection, we estimate the thermal relaxation 
time of a dust particle and then compare it to the rota- 
tion period induced by gas turbulence or by photophorc- 
sis itself. We use a simple model of a cylindrical dust 
particle described in Appendix A. Consider a cylindrical 
dust particle with radius a and height 2a (see Figure IH]). 
The stellar radiation flux / irradiates the front surface. 
We suppose that, at the beginning, the temperature in- 
side the particle is homogeneous. This initial (or average) 
temperature Td is given by equation (|A16[) . assuming 
that the incident flux / balances with thermal radiation 
from the front and back surfaces. (Radiation from the 
side surface is ignored for simplicity.) The temperature 
of the front surface increases to the equilibrium value Tf 
in a thermal relaxation time r t h- During r t h, the incident 
energy at the front surface conducts for a length Sa and 
forms a skin layer of temperature gradient. During rth, 
the temperature at the front surface is approximated as 
the initial value T d , and the energy flux inside the skin 
layer is estimated as Q CO n ~ I — ^sb^I = 1/2. (We 
set e — 1.) When the skin layer depth has grown to 
5a and the temperature of the front surface converges to 
Tf, the energy flux is Q CO n ~ kd(Tf — Td)/5a. Equating 
the above two expressions for Q C on, the thickness of the 
temperature gradient layer is 



Sa 



2k d T d A 



(27) 



where A = (Tf — Td)/Td. Note that the maximum value 
of Sa is 2a. The thermal relaxation time is 



nh 



CdPdSa 2 



(28) 



where Cd is the specific heat capacity of the particle. 

Figure [7] shows the thermal relaxation time r t h for 
10 pm — 10 cm particles in the disk of model A. We set 
the particle bulk density pd = 1 g cm" 3 , the specific heat 
capacity Cd — 10 7 erg g" 1 K -1 , and the thermal con- 
ductivity kd = 10 2 erg s" 1 cm" 1 K _1 . The t emperature 
difference A is calculated by equation (|A17|) . At small 
distances from the star, radiative cooling AeasBTd domi- 
nates in determining A in equation (|A17|) . and Aw 1/4. 



Thus, T t h is independent of the particle size a, and is 
proportional to T^ 6 . In our model Td oc r" 1 / 2 , and thus 
r t h oc 7' 3 . At large distances, on the other hand, internal 
thermal conduction kdTd/a determines A. In this case, 
the skin depth is Sa = a. Hence, r t h is proportional to 
a 2 and independent of r. 

If the rotation period of the particle is smaller than 
T t h, photophoresis is considerably suppressed. Particle 
rotation can be excited by several mechanisms such as 
Brownian motion, gas turbulent motion, collisions with 
other dust particles, and the photophoretic force itself 
can induce rotation if it has a offset from the mass cen- 
ter. Detailed calculation of the rotation period for each 
mechanism is beyond the scope of this paper. Here, we 
make a rough estimate of the rotation speed induced by 
turbulence and by the photophoretic force. 

Suppose that the gas disk has isotropic turbulence 
with the Kolmogorov energy spectrum. We assume that 

the largest eddies have size L 

_ 1/2 
^tur c sj 



a\^h g and velocity 



V = a^c s , where atur is the "a- viscosity" parame- 
ter (Cuzzi et al. 2001). The energy of turbulent mo- 
tion cascades down to smaller eddies, and finally dis- 
sipates by the molecular viscosity. From dimensional 
analysis, the energy dissipation rate per unit mass is 
e ~ V 3 /L ~ atm-Cs^K- The eddy turnover time is faster 
for smaller eddies and that for the smallest eddies is 



Ted 



p g e 



1/2 



-1/2 
*tur 



1/2 



(29) 



where rj V i S = lvxPg/2 is the molecular viscosity (Weiden- 
schilling 1984). The smallest eddies can induce a parti- 
cle rotation period as short as r e d, if the particle is well 
coupled to the turbulent motion of the smallest eddies. 
However, if the particle does not strongly couple to the 
gas, the rotation period is longer than r 0( j. In Figure [3 
the turnover time, r e d, of the smallest eddies is plotted 
by the dashed line for a t ur = 10" 2 . The particle rotation 
period induced by gas turbulence is expected to be larger 
than the dashed line, and therefore is much longer than 
the thermal relaxation time. 

We next consider rotation induced by the torque ex- 
erted by photophoresis. Consider again a cylindrical dust 
particle and suppose that the photophoretic force is ex- 
erted on a point at a distance b from the cylinder axis. 
The torque is K — bF p h. The value of b is unknown, and 
we treat & as a free parameter. The principal moment of 
inertia of the cylinder (for the axis perpendicular to the 
cylinder axis) is I xx = 7ira 5 pd/6. The time needed for a 
initially stationary particle to rotate 180 degree is 



271-4; 
K 



1/2 (7^p d a*\ 1/2 
V 3&F ph 



(30) 



In Figure [51 the rotation time r rot is plotted for assumed 
off-centers b = a and b = 0.01a. If the photophoretic 
force has a large off-center (b = a), the rotation time can 
be shorter than the thermal relaxation time at 0.1 — 2 
AU for 100 p,m particles and at 4 — 40 AU for 1 cm 
particles. In such regions, photophoresis is probably sig- 
nificantly suppressed. In order for photophoresis to work 
effectively in the whole disk, the off-center of the pho- 
tophoretic force has to be as small as 6 = 0.01a. 
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Fig. 8. — Times taken for a particle to rotate 180 degree by 
photophoresis are plotted by the dashed lines. The photophoretic 
force is exerted on a point at a distance b = a or b = 0.01a from 
the cylinder axis. The solid line shows the thermal relaxation time, 
(a) For 100 jim particles, (b) For 1 cm particles. 

5.2. Characteristic Structure Due to Photophoresis 

The dust disk structure that photophoresis makes has 
three zones: the outer disk, the inner disk, and the inner 
hole. The outer disk is composed of small (a < 1 mm 
or si 0.01) particles, and their dynamics is controlled 
mainly by radiation pressure. The inner disk is composed 
of large (a > 100 fim or (3 < 0.01) particles accumulating 
there due to photophoresis. The boundary between the 
outer and inner regions is at 10 — 100 AU, depending 
on the gas density profile. The inner hole opens inside 
~ 0.1 AU, where the Knudsen number Kn is smaller than 
unity for particles that photophoresis effectively works on 
(a = 100 fxm — 10 cm). 

The structure formed by photophoresis should be com- 
pared to the structure of gas-free disks. Modeling of gas- 
free dust disks also shows zonal structure that consists 
of the outer extended disk, the inner disk, and the inner 
hole (Thebault & Augereau 2005; Wyatt 2006; Strubbe 
& Chiang 2006; Krivov et al. 2006). In their models, 
a planetesimal belt is assumed and dust particles are 
continuously produced by planetesimal collisions. The 
inner disk is composed of large particles whose orbits are 
hardly affected by radiation pressure and are nearly cir- 
cular. Thus, the location and width of the inner disk are 
basically similar to those of the planetesimal belt. Out- 



side the planetesimal belt, the outer disk extends and 
it is composed of small particles whose /3-value is large. 
Their orbits are strongly influenced by radiation pres- 
sure and are excited to highly eccentric orbits. In the 
model of Thebault & Augereau (2005), the particles of 
the outer disk have (3 > 0.05, which is relatively larger 
than j3 of our models. In our models, particles' j3 of the 
outer disk can be as small as 0.01. Therefore, infrared to 
millimeter radio observations of outer disks and determi- 
nation of the particle size will provide key information 
for determining which model is plausible. 

Our models of photophoresis also predict that an inner 
hole opens in the dust disk, but the gas still fills the 
hole. This is a significant feature, but it is difficult to 
test this with the present observational techniques. We 
have to observe the gas, but the amount of the gas of 
optically thin disks is probably quite small, and most of 
Vega-type stars do not have gas in an amount exceeding 
the current detection limit. Further, the radius of the 
inner hole is of the order of 0.1 AU, and the temperature 
there is close to the sublimation temperature of the dust. 
(The temperature exceeds 1500 K inside 0.03 AU of our 
model disks.) Therefore, we need a careful observation 
that can distinguish between the inner holes made by 
photophoresis and by dust sublimation. 

5.3. Summary 

Dust accumulation by photophoresis is studied. We 
use formulae of the photophoretic force that are appli- 
cable for the free molecular regime and for the slip-flow 
regime. The main results are as follows. 

1. Particle accumulation occurs at a point where the 
outward acceleration on the gas by the pressure gradi- 
ent equals to the outward acceleration on the particle by 
radiation pressure and photophoresis. 

2. Photophoresis makes an inner disk composed of 
relatively large particles (a = 100 fim — 10 cm). The 
inner disk extends between 0.1 AU and 10 — 100 AU, 
surrounded by the outer disk composed of small particles 
(a < 1 mm). 

3. An inner hole opens inside ~ 0.1 AU. The inner 
hole radius is determined by the condition Kn = 1 for 
the maximum size particles that photophoresis effectively 
works on (a = 100 fim — 10 cm). 

Photophoresis works effectively only when the disk is 
optically thin. Most of small (< 10 /im) dust grains must 
be removed from the disk such that their column density 
to the star becomes smaller than 10~ 3 g cm~ 2 . For ex- 
ample, at 1 AU, the dust density must be smaller than 
10~ 16 g cm" 3 , i.e., 10~ 5 times smaller than the value of 
the minimum mass solar nebula model (Hayashi et al. 
1985). The gas disk also must be optically thin. Fig- 
ure [5] shows that the gas surface density must be smaller 
than 10 2 — 10 3 g cm~ 2 . At 0.1 AU, where the typical ray 
from the star enters the disk (for the dust particles at 
r ^> 0.1 AU), this value is ~ 10 -2 times smaller than the 
value of the minimum mass solar nebula model (Hayashi 
et al. 1985). Even in such a tenuous gas disk, the pho- 
tophoretic force is strong enough to change the dust disk 
structure inside a few AU, as shown by model A (Figure 
Wp)- If the gas density is more tenuous than model A, 
the region where photophoresis has a substantial effect 
shrinks toward the central star. In a gas disk with 0.1 
times smaller gas density than model A, photophoresis 
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Fig. 9. — A cylindrical dust particle of radius a and height 2a. Incident stellar flux / is parallel to the cylinder axis (2-axis) and irradiates 
the front surface. Radiation fluxes from the front and back surfaces are ecrg^T^ and eag^T^ , respectively. The fluxes of heat conduction 
to the gas are qf and qi,, and the heat flux inside the particle is Cj/con- 

on millimeter sized particles works effectively only inside referees for their constructive comments. This work was 
1 AU, and in a disk with 10~ 2 times smaller gas density, supported by the 21st Century COE Program, "The Ori- 
the effective region shrinks to 0.3 AU. gin and Evolution of Planetary Systems," of the Ministry 

of Education, Culture, Sports, Science, and Technology 
„ T ,, , T -j at i i • a t i • ■ t ■ Ti i (MEXT) of Japan, and also by the Grant-in- Aid for Sci- 

We thank Ingrid Mann Tadashi Mukai, Yoichi Itoh, ^ tific ^search, No. 17740107 and No. 17039009, of 
and Yoshitsugu JNakagawa for useful discussions. We are a tut 1 • A c 1 A r< r, j w a 

j. £ i j. rni i • t? •• r j. ■ jiirij- the ME X 1 . OK is grateful to Gerhard Wurm and ac- 

gratetul to lakasm Isuu tor extensive and nelprul dis- , , , r ^ , , ^ , 

° . t r i • 11 -j- i • -ii knowledges the support of the Deutsche 1 orschungsge- 

cussions and for kindly providing his numerical code to meinschaft 
calculate the gas opacity. We are grateful to anonymous 

APPENDIX 

A SIMPLE ESTIMATE OF THE PHOTOPHORETIC FORCE 

In this Appendix, we estimate the magnitude of the photophorctic force using a simple calculation. 

We consider a cylindrical dust particle with radius a and height 2a (see Figure [5]). The particle is surrounded by 
the gas of temperature T g . In general, T g differs from the average temperature of the particle T^. In deriving the 
photophorctic force below, we assume T g = Td, and consider T g ^ Td cases later. The incident radiation flux, /, is 
parallel to the cylinder axis (z-axis) and irradiates the front surface. The temperature at the front surface increases 
to Tf and it radiates energy flux easBTj, where e is the emissivity. Some of the incident energy is conducted to the 
gas with a rate qf per unit area, and the residual energy flows to the back surface with a rate Q C on per unit area. For 
simplicity, we ignore the energy transfer through the side surface. The temperature of the back surface is The 
back surface radiates energy flux eassT^, and conducts heat to the gas with a rate qb ■ In an equilibrium state, the 
temperature inside the particle varies linearly along the cylinder axis, and thus its values at the front and back surfaces 
are written as 

7>=T d (l + A), (Al) 
T b = T d (l-A), (A2) 

respectively. The energy balance reads 

I - sasBTf — qf = Qcon = £tT SB T b 4 + qb . (A3) 

We calculate each term of equation (|A3|) . First we consider free molecular photophoresis (Kn ^1). The radiative 
fluxes from the front and back surfaces are, assuming A <C 1, 

sa SB Tf = £a SB T d 4 (l + 4A) , (A4) 
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ea SB T* = ea SB T$(l - 4A) , (A5) 

respectively. The thermal conduction to the gas at the front surface is divided as 

IS =Qj ~9f , (A6) 

where qj is the energy loss from the particle that is taken by the gas molecules ejected from the front surface, and q^j 
is the energy gain that is given by the adsorbed molecules. In the free molecular approximation (Kn ^> 1), the velocity 
distribution of the adsorbed molecules is written as af + (v), where < a < 1 is the accommodation coefficient, f + (v) 
is the Maxwellian velocity distribution of the temperature T g and number density n g , 

and m g is the mass of a molecule. We define f + (v) to have non-zero values even for particles going away from the 
surface (v z > 0). The energy gain is 



f 1 1 

qj = a -m g v z \v\ 2 f + dv = -aPv T 

Jv z <0 2 2 



(A8) 



where v z is the velocity component normal to the front surface, P — npgV^/8 is the gas pressure, and vt = 
yJ&kBTg/itrrig is the mean thermal speed of gas molecules. The molecules ejected from the front surface are assumed 
to have a Maxwellian velocity distribution, af~(v), of the temperature Tf. The energy loss is 



r i n 

qj =a ^m g v z \v\ 2 f-dv = T^ap](vj) 3 , 



(A9) 



where — ^/SksTf /irm g and the "density" of ejected molecules (for a = 1) is 

Pj = m g f-dv . (A10) 



all V 



From mass conservation between the adsorbed and ejected molecules {oipj v j = apgVr), 



T, 



1/2 



Then, we obtain 



Pj = 1 ) 0, ■ (All) 

<7 ~ ■ (A12) 

If the gas temperature and the average dust temperature are same (T g — Td), the thermal conduction to the gas is 

q f = ^aPv T A . (A13) 

Similarly, the thermal conduction to the gas from the back surface is 

q b = -^aPv T A . (A14) 

The thermal conduction inside the particle is 

Qcon = . (A15) 

a 

Collecting all the above expressions of the terms in equation (|A3[) . the particle temperature is solved as 

r *=(^) ' (A16) 

and 

A = - (A17) 

2[k d T d /a + 4e<j SB Tl + aPv T /2} ' { ' 

The force exerted on the front surface by the ejected molecules is 



r 

Ff = ira a / m g v z f~dv = 

Jv z >Q 



ira 2 aP (T^ 1/2 



(A18) 
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and on the back surface, 

Assuming T s = T<j and A< 1, the photophoretic force is 



Tra 2 aIP 

ph ' f = ' 6 = IwTwfT^ ' (A20) 

which can be compared to equation for the spherical particle case. 

We consider how the photophoretic force changes from equation (|A20[) when the dust temperature Td differs from the 
gas temperature T g . For simplicity, we consider the cases where the thermal conduction to the gas can be neglected. 
With qj = qb = in equation (|A3[) , the photophoretic force becomes 

JW __f^_^V /, 1 (A21) 



4[k d T d /a + AsctsbT^ \T g , 

which differs from the T d = T g case by a factor (Tj/Tg) 1 / 2 . We consider further the cases where the temperature 
variation inside the dust particle is large (A ~ 1). Such a situation arises if easBT^ S> kdTd/a, aPvr- In this case, 
Tf = (I/ea S B) 1/4 and T b =0, and thus the photophoretic force is F ph)f = F f = na 2 aP(T f /Tg) 1 ' 2 /2. The photophoretic 
force calculated by equation (|A21|) . which assumes A <S 1, r educe s to i^ p h,t = Ff/{A ■ 2 1 / 8 ), when the term kdTd/a is 
neglected. Thus, in the extreme case where A ~ 1, equation (|A21[) underestimates the photophoretic force by a factor 
4-2 1 /s. 

In the slip-flow regime (Kn <C 1), thermal conduction to the gas is estimated as qj ~ qb ~ k g TdA/a, where k g is the 
thermal conductivity of the gas. The temperature difference is 

A ~ - (A22) 

[{kd + k g )T d /a + 4ea S BT^} ' V ' 

Because of the temperature variation along the dust surface, a slip flow of the gas arises. The velocity of the thermal 
slip is 

v s - C s « r AKn , (A23) 

where C s ~ 1 is the thermal slip coefficient (Lifshitz & Pitaevskii 1981). The particle, which moves in the gas with a 
velocity v s , experiences Stokes gas drag, 

-Fstokcs ~ aihisVs (A24) 

where ?7vis = lvrPg/2 = VTm g / (2^/2a mo i) is the molecular viscosity. (-Fstokcs = &irarj w - ls v s if the particle is spherical.) 
In the equilibrium state, the photophoretic force balances with the Stokes drag force, and is estimated as 

C s a 2 k B I A , 4aeo- SB T| . kA' 1 

Fph.s : — 1 + : + 7^ Kn , (A25) 

v mo \kd \ k d k d J 

which can be compared to equation ([5])- In the above derivation, we neglect the temperature jump at the surface and 
the velocity slip due to imperfect sticking of the gas to the surface (Ct = C m = 0). 

SLIP-FLOW REGIME PHOTOPHORESIS 

In this Appendix, we derive the expression for the photophoretic force in the slip-flow regime, equation (|5|), taking 
into account the radiative cooling of the particle. We follow the derivation described in Mackowski (1989). (The 
assumption is as follows. The mean temperatures of the gas and the dust, T g and T^, are the same. The temperature 
variation from the mean value is small. The Reynolds number is much less than unity and thus gas drag obeys the 
Stokes law. The gas temperature is determined only by the heat conduction from the dust and inside the gas.) Here, 
we do not rewrite Mackowski's derivation again, but we describe what is modified when radiative cooling is taken into 
account. In this Appendix, the numbering of equations with "M" means that they refer to equations of Mackowski 
(1989). For the meaning of the symbols see Mackowski (1989). However, for some variables, we have used different 
symbols in the main text, and we keep the same notation in the Appendix. When we use a different symbol from 
Mackowski (1989), the corresponding symbol in Mackowski (1989) is given in brackets. 

The basic equations (M11)-(M14) and the boundary conditions (M15)-(M20) are not changed except that equation 
(M16) is modified as 



k g ^ + ea S BTd 4 = ~k d ^ , (Bl) 
or or 



where Td(T s ) is the temperature inside the dust particle, T g (T g ) is the gas temperature, and fcd(fc s ) is the thermal 
conductivity of the dust. Then, equation (M22) becomes 

D = nG n {\) - G' n {\) 

" (Ae<j SB T^a/k d + n)[l + (n + l)C t l/a} + (n + i)k g /kd 1 1 ' 
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where Id(Xo) is the mean temperature of the dust, Ct{c t ) is the coefficient of the jump condition, and we correct the 
typo in the sign of the original equation. Equation (M28) is not modified when it is expressed using D\ as 

= 3(1 + 2C m l/a) _ CsT/yisDi 
Cl 4(l + 3C m //a) 2V oPg a{l + 3C m l/a) ' [ ' 

where r] v i s (r]) is the molecular viscosity, and C m (c m ) and C s (c s ) are the coefficients of the jump conditions. Finally 
equation (M29) becomes 



= AirC^aLh 1 

p fc d T d p g (l + 3C m Z/o)[(4e(7-sBr|o/A (J + l)(H-2C t Z/o) + 2fc fl /fcJ ' 



where is the incident flux of the starlight. Substitution of the expression for rfc is — \f2lp g kBTd/ '{^(y mo \) and 

taking the absolute value give equation ([5]). 



REFERENCES 



Alexander, R. D., Clarke, C. J., & Pringle, J. E. 2006, MNRAS, 
369, 216 

Alexander, R. D., Clarke, C. J., & Pringle, J. E. 2006, MNRAS, 
369, 229 

Alexander, D. R., & Ferguson, J. W. 1994, ApJ, 437, 879 
Ardila, D. R., et al. 2005, ApJ, 627, 986 

Augereau, J. C, Lagrange, A. M., Mouillet, D., & Menard, P. A. 

1999a, A&A, 350, L51 
Beresnev, S., Chernyak, V., & Fomyagin, G. 1993, Phys. Fluids A, 

5, 2043 

Borysow, A., Jorgensen, U. C, & Zheng, C. 1997, A&A, 324, 185 
Bottke, W. F., Jr., Vokrouhlicky, D., Rubincam, D. P., & Broz, M., 
in Asteroids III, ed. W. F. Bottke Jr., A. Cellino, P. Paolicchi, 
& R. P. Binzel, (Tucson: University of Arizona Press), 395 
Brownlee, D., et al. 2006, Science, 314, 1711 
Burns, J. A., Lamy, P. L. & Soter, S. 1979, Icarus, 40, 1 
Chapman, S., & Cowling T. G. 1970, The Mathematical Theory of 

Non-Uniform Gases (Cambridge: Cambridge University Press) 
Chen, C. H., & Kamp, I. 2004, ApJ, 602, 985 
Clampin, M., et al. 2003, AJ, 126, 385 

Clarke, C. J., Gendrin, A., & Sotomayor, M. 2001, MNRAS, 
328,485 

Coulson, I. M., Dent, W. R. F., & Greaves, J. S. 2004, MNRAS, 
348, L39 

Cuzzi, J. N., Hogan, R. C, Paque, J. M., & Dobrovolskis, A. R. 

2001, ApJ, 546, 496 
Dent, W. R. F., Greaves, J. S., & Coulson, I. M. 2005, MNRAS, 

359, 663 

Dullemond, C. P., & Dominik, C. 2005, A&A, 434, 971 
Ferguson, J. W., Alexander, D. R., Allard, F., Barman, T., 

Bodnarik, J. G., Hauschildt, P. H., Heffner-Wong, A., & 

Tamanai, A. 2005, ApJ, 623, 585 
Fisher, R. S., Telesco, C. M., Pina, R. K., Knacke, R. F., & Wyatt, 

M. C. 2000, ApJ, 532, L141 
Goto, M., Usuda, T., Dullemond, C. P., Henning, Th., Linz, H., 

Stecklum, B., & Suto, H. 2006, ApJ, 652, 758 
Greaves, J. S. 2004, MNRAS, 351, L99 

Greaves, J. S., Mannings, V., & Holland, W. S. 2000, Icarus, 143, 
155 

Hartmann, L., Calvet, N., Gullbring, E., & DAlessio, P. 1998, ApJ, 
495, 385 

Hayashi, C, Nakazawa, K., & Nakagawa, Y. 1985, in Protostars 
and Planets II, ed. D. C. Black & M. S. Matthews (TucsomUniv. 
of Arizona Press), 1100 
Hidy, G. M., & Brock, J. R. 1967, J. Geophys. Res., 72, 455 
Hollenbach, D., Johnstone, D., Lizano, S., & Shu, F. H. 1994, ApJ, 
428, 654 



Krauss, O., & Wurm, G. 2005, ApJ, 630, 1088 

Krauss, O., Wurm, G., Mousis, O., Petit, J. -M., Horner, J., & 

Alibert, Y. 2007, A&A, 462, 977 
Krivov, A. V., Lohne, T., & Sremcevic, M. 2006, A&A, 455, 509 
Lifshitz, E. M., & Pitaevskii, L. P. 1981, Physical Kinetics (Oxford: 

Butterworth-Heinemann) 
Liseau, R., & Artymowicz, P. 1998, A&A, 334, 935 
Mackowski, D. W. 1989, Int. J. Heat Mass Transfer, 32, 843 
Marsh, K. A., Silverstone, M. D., Becklin, E. E., Koerner, D. W., 

Werner, M. W., Weinberger, A. J., & Ressler, M. E. 2002, ApJ, 

573, 425 

Miyake, K., & Nakagawa, Y. 1993, Icarus, 106, 20 
Mouillet, D., Lagrange, A. M., Augereau, J. C, & Menard, F. 2001, 
A&A, 372, L61 

Mousis, O., Petit, J.-M., Wurm, G., Krauss, O., Alibert, Y., & 

Horner, J. 2007, A&A, 466, L9 
Nakagawa, Y., Sekiya, M., & Hayashi, C. 1986, Icarus, 67, 375 
Nomura, H., & Nakagawa, Y. 2006, ApJ, 640, 1099 
Petit, J.-M., Mousis, O., Alibert, Y., & Horner, J. 2006, Lunar 

Planet. Sci., 37, 1558 
Presley, M. A., & Christensen, P. R. 1997, J. Geophys. Res. E, 102, 

6535 

Strubbe, L. E., & Chiang, E. I. 2006, ApJ, 648, 652 

Takeuchi, T., & Artymowicz, P. 2001, ApJ, 557, 990 

Takeuchi, T., Clarke, C. J., & Lin, D. N. C. 2005, ApJ, 627, 286 

Tanaka, H., Himeno, Y., & Ida, S. 2005, ApJ, 625, 414 

Thebault, P., & Augereau, J.-C. 2005, A&A, 437, 141 

Tsuji, T. 1966, PASJ, 18, 127 

Tsuji, T. 1971, PASJ, 23, 553 

Tsuji, T. 2002, ApJ, 575, 264 

Weidenschilling, S. J. 1984, Icarus, 60, 553 

Weidenschilling, S. J., & Cuzzi, J. N. 1993, in Protostars and 

Planets III, ed. E. H. Levy & J. I. Lunine (TucsomUniv. of 

Arizona Press), 1031 
Weinberger, A. J., Becklin, E. E., Schneider, G., Smith, B. A., 

Lowrance, P. J., Silverstone, M. D., Zuckerman, B., & Terrile, 

R. J. 1999, ApJ, 525, L53 
Weinberger, A. J., Rich, R. M., Becklin, E. E., Zuckerman, B., & 

Matthews, K. 2000, ApJ, 544, 937 
Wurm, G., & Krauss, O. 2006, Icarus, 180, 487 
Wyatt, M. C. 2006, ApJ, 639, 1153 

Zuckerman, B., Forveille, T., & Kastner, J. H. 1995, Nature, 373, 
494 



